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ABSTRACT 

The concordance particle creation model - a class of A(t)CDM cosmologies - is studied 
using large scale structure (LSS) formation, with particular attention to the integrated 
Sachs-Wolfe (ISW) effect. The evolution of the gravitational potential and the ampli¬ 
tude of the cross-correlation of the cosmic microwave background (CMB) signal with 
LSS surveys are calculated in detail. We properly include in our analysis the peculiar¬ 
ities involving the baryonic dynamics of the A(t)CDM model which were not included 
in previous works. Although both the A(t)CDM and the standard cosmology are in 
agreement with available data for the CMB-LSS correlation, the former presents a 
slightly higher signal which can be identified with future data. 
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1 INTRODUCTION 


pm -f pm 


-A, 


( 1 ) 


Although the great success of the standard ACDM cosmo¬ 
logical model in describing most observations, we are still 
distant from the full understanding of the cosmic dynamics. 
Apart from the so called small scale problems in the dark 
matter sector, there is also a huge discussion about dark 
energy properties. At the same time, the efficient approach 
given by a cosmological constant A still seems to face some 
challenges, mainly from the theoretical point of view. 

Among the viable alternatives, it has been argued that 
a model with the mechanism of dark matter particle produc¬ 
tion at a constant rate F (implying in a dynamical vacuum 
term A{t)) is also capable to produce a concordance model 
l Alcaniz et al.|2012 l, although with a larger current matter 
fraction QmO ~ 0.4E^ 

If the total energy is given by p = pm + A and the 
pressure of the vacuum contribution sets to pA = —A, a 
constant F is equivalent to take A = 2FF/, where H = a/a 
is the Hubble rat^ Taking today’s values one finds F = 
3 DaoHo/ 2. The matter dynamics corresponds to the balance 
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^ For any component i we define ULi = pi/pcO with p^o = 
3Hq /SttG, where the subscript 0 denotes today’s values. In this 
paper we will take SttG = c = 1. 

^ Strictly speaking, this equivalence is only valid if we neglect the 
conserved baryons in the background equations. In the presence 
of a small baryonic content, taking A oc R does not lead to an 
exactly constant F. 


from which it is not difficult to show that A = —Fp^. The 
background expansion is given by 

nmo + , ( 2 ) 

from which we can derive the matter density evolution 

Pm = 3iFo “f (1 — flmo)^mO^ ^ ] • (3) 

Note that, as in the flat ACDM model, there are here also 
only two free parameters. The concordance of this A(t) cos¬ 
mology with VlmO ~ 0.45 has been verified via many different 


observational data at both the background 

Carneiro et al. 

2006 2008| Pigozzo et al.|20Il 

Velten et al. 

2005 \ and per- 

turbative (|Alcaniz et al.||2012 

Zimdahl et al.||20Il Devi et 


Concerning structure formation, the full CMB spectrum 
has not yet been obtained for this model, but CMB physics 
can also be accessed via ISW studies, i.e., how time-varying 
gravitational potential wells change the temperature of the 
CMB photons as they cross structures (Sachs et al.||1967|. 


The expansion of an Einstein-de Sitter universe compensates 
the clustering of structures, producing no ISW effect, i.e., in 
a matter dominated universe there is no “late time” ISW ef- 
fecl[^ Dark energy modifies the background expansion and 


® The “early time” ISW effect is related to non-insignificant ra¬ 
diation density just after photon decoupling. 
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leads to a net contribution to (AT/T)^^^. Then, it is ex¬ 
pected that the modified background and perturbative ex¬ 
pansion of the A(t)CDM model leaves a distinct imprint on 
the ISW signal. 

Since the ISW is a secondary CMB temperature effect 
its detection occurs only via the cross-correlation with other 


large scale probes like galaxies and quasars surveys (Crit- 
tenden et al.||2005 |. Ref. ( Wang et al.|2010| used the cross¬ 


correlation technique to probe the A(t)CDM model, finding 
an increase of the ISW-galaxy spectrum iCf^) in compari¬ 
son to the ACDM model. Our aim in this paper is to per¬ 
form this analysis taking care with some peculiarities of the 
model. In particular, it is important to differentiate the evo¬ 
lution of the baryonic and dark matter components. At the 
background level, baryons are included in the conserved part 
of §. Concerning perturbations, the observed matter power 
spectrum P{k) = |5(fc)|^ is a measurement of the baryonic 
clustering, which is sourced by the gravitational potential. 
In the ACDM model there is almost no difference between 
the evolution on large scales of such components. The same 
does not happen in the A(t)CDM model, where one observes 
a late-time suppression in the dark matter and total matter 
contrasts owing to dark matter production, a suppression 
not observed in the baryonic contrast. 

Our interest is also justified by a possible tension be¬ 
tween the theoretical ACDM predictions and the observed 
ISW signal. Some analysis have reported a cross-correlation 
signal above that expected in the standard cosmol¬ 


ogy (Giannantonio et al. 2008 Goto et al. 20121 (see (Gi- 
annantonio et al.||2012l for a critical review on this issue). 


Recent studies claim a better statistical concordance though 
the observed signal is still higher than the theoretical one 
( Kovacs et al.|[^13 Ferraro et al.|[20T4 l. Also, the inferred 
stacking of CMB data on the position of superstructures is 


about 5 times larger than the ACDM prediction ( 

Grannett 

et al.|2008| INadathur et al.|20121|Hernandez-Monteagudoet 

al. 2013 

Flender al. 20131. This result has recently been con- 


firmed by the 2015 release of the Planck CMB data (Planck 
|collaboration|[2012[ ). 

In the next section we develop the background expan¬ 
sion which will be used in this paper. In sectionj^we explore 
the scalar perturbations of the model. The connection with 
the ISW effect is made via a direct calculation of the evo¬ 


lution of the gravitational potential (3.21 and the ISW-LSS 


cross-correlation (13.31). We conclude in the final section. 


2 A(T) BACKGROUND DYNAMICS 

The brief description of the A(t)CDM dynamics given in 
equations ijs has been widely derived in the literature (see 
(Alcaniz et al. 2012) and references therein). Note that as 


in the standard case the important quantity is the total 
matter density DmO, the sum of the dark matter fldmo and 
baryons The latter is constrained by nucleosynthesis 
and is well approximated by = 0.05. We assume here¬ 
after this value. The baryonic sector is conserved and there¬ 
fore we have the continuity equations 


pdm pdm — A 

and 

pb + 3Hpb = 0. 


( 4 ) 


( 5 ) 


Adding these equations we obtain the total matter conserva¬ 
tion equation (Q. With A = 2rjT, the dynamics presented 
in section follows. That is, the background expansion is 
obtained after solving the system given by the continuity 
equation 

aHp„, + 3Hpm = Ppm (6) 

and the Friedmann equation 

2aHH' -b = 0, (7) 

where the prime means derivative with respect to the scale 
factor. 


3 PERTURBATIVE DYNAMICS 
3.1 Growth functions 

When cosmological perturbation theory is used to study the 
large scale structure of the universe we are most interested 
in the evolution of scalar quantities like the density con¬ 
trast and the gravitational potential. The A(t)CDM mod¬ 
els have the particularity of suppressing the growth of dark 
matter contrast 5dm- However, the observed structures re¬ 
flect the behavior of the baryonic contrast, 5b = 5ph/pb, 
instead. In the standard cosmology this quantity coincides 
with the dark matter contrast, and the total matter con¬ 
trast 5m denotes both behaviors. For the A(t)CDM model, 
however, we should make a distinction between the different 
components. Using the fact that the vacuum fluctuations are 
negligible ( Zimdahl et al.|2011 1 and choosing the comoving 
synchronous gauge we obtain the set of equation^ 


.+V5m = ^, 


( 8 ) 


5b = 


2 ’ 


h -b 2Hh = pm5r, 


( 9 ) 

( 10 ) 


where h is the trace of the spatial metric perturbation. Gom- 
bining these equations and changing to the derivative with 
respect to the scale factor we find 


+ aH{ZH -b aH' -b r)C + 


2Y^H5m — i^pm5m 


( 11 ) 


and 


+ aH{SH + aH')5'b = -pm5r, 


( 12 ) 


In Fig. we plot the evolution of the density contrasts 


as calculated in (111 and (121. For the AGDM model (black 


curves) we set 

„o = 0.23, the best-fit of LSS observa- 

tions ( 

Cole et al. 

2005 

Percival et al. 

2007f 1 We show 


the A(t)GDM evolutions for Sb and 5m in the solid red 
and dashed red curves, respectively. We have used initial 
conditions 5b{a = 0.001) = 5m{a, = 0.001) = 10“® and 
5'b(a — 0.001) = 5m(a- = 0.001) = 0. As expected, 5b and 


A different approach to the baryonic sector can be found in 
ijZimdahl et al.|2014|. 

® Using, instead, the ACDM concordance value fimO = 0.27 does 
not change significantly our final results. 
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a 

Figure 1. Evolution of the linear density contrasts. The black 
curves correspond to the baryonic (solid-black curve) growth and 
total matter (dashed-black curve) growth for the ACDM model 
with flrnO = 0.23. Related plots for the A(t)CDM with QrnO = 
0.45 are in solid red (baryons) and dashed red (total matter), 
respectively. 


Sm have the same evolution in the standard case. On the 
other hand, the plot shows a late-time growth suppression 
of Sm in the A(t)CDM model, a consequence of dark mat¬ 
ter creation. The same suppression is not observed in the 
baryonic contrast though its amplitude achieves a plateau 
on late times which is slightly below the standard scenario. 


3.2 Evolution of the gravitational potential 


The evolution of the gravitational potential <0 determines 
the integrated Sachs-Wolfe contribution to the total CMB 
spectrum via the formula 


^rpISW 

~T~ 




(13) 


where the subscript “d” denotes the decoupling time, and 
a = 1 at present. Writing down the Poisson equation in the 
comoving synchronous gauge we find 
2 

k^^ = -^PmSm, (14) 


where k is the comoving wave-number. This expression al¬ 
lows us to calculate the ISW, which can be computed via 
the line of sight integration ( |13[ ). It is useful to perform it in 
terms of the comoving distance 


X{a) 


e da 

Ja o?H{ay 


(15) 


The ISW then reads 

= ^ dQm{a), 

-T- 

where we have defined 


(16) 


Qm (Q)) 

with 

D^ia) 


a^Pm{a)Dy{a), 

Sm (a) 

Sm(a = 1)' 


(17) 

(18) 


The baryonic and total matter density contrasts can be 
directly calculated from the solution of Eqs. (Ill and (12|. 



z 


Figure 2. Evolution of the gravitational potential for the ACDM 
model with QmO = 0.23 (black) and for the A(t)CDM with QmO = 
0.45 (red-dashed). 


Using Eq. ( |14[ ), we compare in Fig. I^the predictions for the 
ACDM and A(t)CDM gravitational potentials as functions 
of the redshift z. We access the observational predictions for 
the ISW effect via the cross-correlation of CMB maps and 
LSS surveys, to be done in next section. 


3.3 Correlating CMB maps and galaxy surveys 


In order to compute the cross-correlation between CMB and 
LSS we also need to describe the evolution of the observed 
galaxy contrast Sg on the line of sight. This quantity depends 
on the survey design and is obtained from 

^9=J b{z)-^H{a)Sb{a)dx, (19) 


where Si, is the baryonic linear contrast, and b{z) is the stan¬ 
dard bias between Sg and 

The redshift distribution of the observed galaxy sam¬ 
ple is a model-independent quantity. Each survey has its 
own dN/dz histogram function. We use in this work data 
from the NRAO VLA Sky Survey (NVSSQ and the Wide- 
field Infrared Survey (WISE)[^ catalogues. Both surveys are 
widely used for cross-correlation studies. The NVSS covers 
the entire north sky of —40 deg declination in one band. The 
produced catalogue of discrete sources contains more than 


one million objects. Full details appear in Ref. (Condon et 
al.|1998 1. WISE has an entire sky scanning strategy in four 


frequency bands. It has an average redshift 0.3, but reach¬ 
ing up to a ~ 1, collecting about 500 million sources which 
are in general not restricted to point-like objects such as, 
for example, stars and unresolved galaxies. The data used 
in this work is taken from the analysis recently performed 
in Ref. (( Ferraro et al.|2014 |) which uses a larger sample by 
applying less-conservative cuts to the dataset in comparison 
to previous works ( Goto et al.|2012 Kovacs et al.|20l3 |. 

For the NVSS catalogue one has 


® In the ACDM model is irrespective to use 5^ or Sm in | |19[ |, but 
the same is not true for the A(t) case, as can be seen in Fig. 

http://www.cv.nrao.edu/nvss/ 

* http://wise.ssl.berkeley.edu/ 
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( 20 ) 


o{z)——=beff - n -exp 

^ ' dz " z:+^r{a) ^ 

where the values b^ff = 1.98, 2 * = 0.79 and a = 1.18 were 


az 

2 * 


htted in (Ho et al. 20081. For the WISE catalogue dN/dz 
can be obtained numerically from Ref. | Yan et al.|2010 1 and, 
following Ref. I Ferraro et al.|2014 1, we adopt a constant bias 
6 ( 2 ) = 1.41. We will limit our results to constant bias models 
since the use of time-dependent 6 ( 2 ) bias will not change our 
main conclusion. 


Combining (16 1 and (191, the multipole coefficients for 
the cross-correlation spectrum are given by 


CtAI) = 


WT{a)WAa) Pik = l/x) 


da, 


a‘^H{a) P 

where we have defined the weight functions 
,dQ{a) 


Wria) = a^H{a)- 


da 

dN 


WAa) = H{a)b{z)'^Dt{a), 

with 


Dt{a) = 


6b{a) 


( 21 ) 

( 22 ) 

(23) 

(24) 


5b{a = 1)' 

We have also defined the crossed power spectrum 

P{k) = (25) 

di,(a = 1) 

where Pb{k) is the observed baryonic spectrum, given by 


Pbik) ^ Pok"‘TAk), 


(26) 


where Po is a normalisation constant determined from ob¬ 
servations ( Alcaniz et al.|2012[ ). This normalisation constant 
is not the same as in the ACDM model as can be seen by an 
inspection of figures [^and It can also be determined by 
writing Po( 55 (l) = Po^(l), where the barred quantities are 
the standard ones ( |Devi et al.||2015[). For T {k) we will use 
the BBKS transfer function ( Barden et al.|[r986| Martin et 
^|2000 1 which, for fli,o flmo, can be approximated by 


T(x = k/keq) = 


ln[l + 0.171a;] 
(0.171a:) 


(27) 


[1 + 0.284X -b (l.lSa;)^ -b (0.399a;)^ -b (0.490®)'*] . 

Here, k^q is the comoving Hubble radius at the time of 
matter-radiation equality. Fr om ([3| ) it is easy to show that 
it is given by (Alcaniz et al. 2012 Devi et al. 20151 k^q = 
0.073h^flmo Mpc“ , where, here, h = 770/(100 km/s-Mpc) 
and we have set the present radiation density as flijo = 
4.15 X We will adopt h = 0.7 and Us = 1. 

The resulting cross-correlation spectrum is show in Figs. 
[^and|^ For the NVSS data (Fig. a we follow the presenta¬ 


tion of data used in (Wang et al. 


20101 , where the original 


data points of Ref. ( Ho et al.||2008 1 located at Z < 10 and 
I > 70 have been removed due to their high dispersion. The 
WISE data presented in Fig.is taken from (Ferraro et al. 
2014). In both cases, the C; ® spectrum for the A(t)CDM 


model presents a slightly larger signal as compared to the 
standard model. This seems to be a virtue of the particle cre¬ 


ation model since we have learned from Refs. (Giannantonio 


et al. 

2012 

Kovacs et al.|2013| [Ferraro et al.|2014| Grannett 

et al. 

2008 

Nadathur et al.|2012||Hernandez-Monteagudoet 



Figure 3. CMB-LSS cross-correlation spectrum using the NVSS 
data taken from et al.|2008| l (see text for details). The black 
line represents the ACDM with QmO = 0.23. The red-dashed line 
is our inferred spectrum for the A(t)CDM model using QrnO = 
0.45. 



Figure 4. The same as Fig. for the WISE data taken from 
([Ferraro et al.|2014[). 


that models with higher 
power are desirable. Nevertheless, given the large uncertain¬ 
ties in determining the observed values, both models 
remain compatible with data. 


al. 


2013 


Flender al. 


2013 


4 CONCLUSIONS 

The concordance particle creation model is a viable alterna¬ 
tive to the standard ACDM paradigm, with the same free 
parameters, namely QmO and Hq. Since particle creation at 
a constant rate leads to a different background and pertur¬ 
bative dynamics, it is important to investigate specific sig¬ 
natures of this model in comparison to the standard cosmol¬ 
ogy. We provided in this work a direct computation of the 
evolution of perturbed scalar quantities like the matter and 
baryonic density contrasts and the gravitational potential 
for the A(t)CDM model. There is a clear scenario emerging 
in such a cosmology (see Fig. [^, in which the dark matter 
contrast is highly suppressed at late times while the bary¬ 
onic contrast maintains a constant value, as in the standard 
case, but with a slightly smaller amplitude. This dynamics 
does indeed lead to a consistent description of the structure 
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formation data (Alcaniz et al. 2012 Zimdahl et al. 2011 
Devi et al.|2015 I. Our main goal here was to provide a com¬ 


prehensive analysis of the ISW effect in such concordance 
particle creation cosmology. With the results for the pertur¬ 
bative dynamics calculated In section II we computed the 
CMB-LSS cross correlation spectrum. Our care in calculat¬ 
ing the evolution of both baryonic and total matter contrasts 
in detail leads to a full and safe analysis of the Cj" 
trum. 


Tg 


spec- 


Current efforts | 

Giannantonio et aL|201 

2 Kovacs et al. 

20131 

Ferraro et al. 

2014| Grannett et al.| S 

0081 INadathur 

et al.l 

20121 Hernanc 

ez-Monteagudoet al.| 20131 IFlender al. 

2013) in obtaining the observed C/® spectrum are sending 


a clear message: the CMB-LSS signal seems to slightly ex¬ 
ceed the ACDM prediction. Therefore, it is timely to check 
the A(t)CDM outcomes for this cosmological observable. 
For both the NVSS (Fig. and WISE (Fig. data, the 
A(f)CDM model leads to a desirable excess of power in the 
spectrum of order 10-20%. This is in fact a small in¬ 
crease that potentially cannot be sensitive to ISW probes. 
It is worth noting that the available data is insufficient for 
distinguishing cosmological models with a reliable statistical 
confidence. This happens mainly because of the cosmic vari¬ 
ance limits imposed to large scale CMB analysis, that leads 
to a very low signal-to-noise ratio. Here we are restricted to 
a qualitative comparison between the A{t) and standard cos¬ 
mologies. Future data can in principle improve the accuracy 
of the CMB-LSS cross-correlation technique and therefore 
specific features of different cosmological models concerning 
the ISW effect could be compared in more detail. Although 
pure ISW detection techniques do not seem to be a pow¬ 
erful tool in discriminating cosmological models, they can 
eventually complement other cosmological probes. 
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